Quantum Trajectories for Brownian Motion 
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We present the stochastic Schrodinger equation for the dynamics of a quantum particle coupled to a 
high temperature environment and apply it to the dynamics of a driven, damped, nonlinear quantum 
oscillator. Apart from an initial slip on the environmental memory time scale, in the mean, our 
result recovers the solution of the known non-Lindblad quantum Brownian motion master equation. 
A remarkable feature of our powerful stochastic approach is its localization property: individual 
quantum trajectories remain localized wave packets for all times, even for the classically chaotic 
system considered here, the localization being stronger as Ti — > 0. 
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The understanding of the dynamics of open or dissi- 
pative quantum systems is of fundamental importance 
both from a practical and conceptual point of view. The 
archetype of such a system is the standard quantum 
Brownian motion model [Q which describes a particle 
with Hamiltonian H(q,p), coupled to an environment of 
harmonic oscillators (q\,p\) via its position q, such that 
the total Hamiltonian of system and environment reads 



H to t(q,P, Q\,P\) 
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Up to now, in order to determine the time dependent 
dynamics of the open 'system', the standard procedure 
was the derivation of a master equation for the reduced 
density operator, which, for the high temperature case 
considered below, is widely accepted to read 



i[H,pt] 
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where 7 is the damping rate. This master equation is a 
Markov master equation not, however, of Lindblad form 
H and indeed it turns out that it may violate the posi- 
tivity of p t on very short time scales, which has led to an 
ongoing debate about its range of applicability ||. We 
will briefly address this issue later on in this Letter. 

Our new approach to quantum Brownian motion is 
very different and circumvents the derivation of a master 
equation for p t altogether. Instead, we use a stochastic 
Schrodinger equation, derived straight from the micro- 
scopic model ([!]), for pure states ijjt(z) (quantum trajec- 
tories). Our construction recovers the reduced density 
operator as the ensemble mean M[. . .] over many of these 
quantum trajectories, in principle without any approxi- 
mation: 



The mean M[. . .] is taken over the process z t which drives 
the stochastic Schrodinger equation. We are thus able to 
determine p t in a Monte-Carlo sense without an explicit 
master equation for its time evolution. 

Quantum trajectory methods have been used exten- 
sively in recent years, mainly in the quantum optics com- 
munity, due to their numerical efficiency, their intimate 
connection to (continuous) measurement, and their il- 
lustrative power helping to gain physical insight. The 
master equations encountered in quantum optics are of 
standard Lindblad type, for which Markov quantum tra- 
jectory methods are known for some time now: there are 
jump processes Q] and diffusive processes [[| recovering 
the reduced density operator. Despite being maybe the 
best known of all master equations, the Quantum Brow- 
nian motion master equation (||), being not of Lindblad 
form, has so far been excluded from a treatment with 
these powerful methods. 

Only recently the authors managed to extend the quan- 
tum trajectory concept to non-Markovian situations ||, 
more precisely, we were able to determine a stochastic 
Schrodinger equation for the dynamics of a quantum sys- 
tem coupled to a bath of harmonic oscillators as in (0), 
without using the concept of a master equation for p t . 
An alternative approach to non-Markovian quantum tra- 
jectories, more emphasizing the continuous measurement 
point of view, has now also been established 0. 

In its linear version ||, our non-Markovian quantum 
state diffusion (QSD) stochastic Schrodinger equation for 
the quantum Brownian motion model ([]]) takes the form 



hipt(z) = -iH'^t(z) + qzt^t(z) 
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pt = M[\ilH(z))(Mz) 



(3) 



where we assumed a factorized total initial density oper- 
ator p tot = |^o)(V ; o| ® Pt with a pure system state |^o) 
and an environmental thermal density operator pT. The 



1 



influence of the environment on the system is encoded 
in the bath correlation function a(t,s) = (F(t)F(s)) qt 
where F(t) = ^a^a^a(^) is the quantum force in (|l]) 
and zt is thus a complex Gaussian stochastic c-number 
force with correlation M[z*z s ] = a(t,s). In the usual 
high temperature limit kT ^> ft A ^> hcu,h'y, where A is 
an environmental cutoff frequency and a;, 7 are the typi- 
cal system frequency and damping rate, respectively, one 
finds [| 



a(t, s) = 2mjkTA(t - s) + ihmjA(t - s), 



(5) 



where A(t) is a delta-like function decaying on the 
environmental 'memory' time scale A -1 (here we use 
A(t) = f e" A l*l). In (|), the Hamiltonian H' = H(q,p) + 
^mjAq 2 contains an additional potential term that turns 
out to be counterbalanced by a similar term arising from 
the memory integral. 

Eq. (ji)) is exact, i.e. it provides a quantum trajec- 
tory method for Brownian motion for any temperature 
and any distribution of environmental oscillators in the 
model (|l|), i.e. for any a(t, s). In order to compute num- 
bers, however, we have to express the functional deriva- 
tive under the memory integral in (^) in terms of elemen- 
tary operators. In the high temperature limit considered 
here, we simply need to expand in terms of the time delay 
(t-s) 
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where the dots denote terms of the order (t — s) 2 and 
higher, leading to corrections of the order uj/A, 7/A and 
can therefore be neglected (see (9) for a general theory of 
such 'post-Markov' open systems). With (||), the mem- 
ory integral in (^) takes the form 
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where we introduce time dependent coefficients go(t) — 
\ f*dsa(t,s) and 01 (t) = ^ f*ds(t-s)a(t,s). The imag- 
inary part of go (t) will be compensated by the additional 
potential term in H' . The imaginary part of gi(t) gives 
rise to damping. The real part of go(t) describes diffusion 
and as the real part of gi(t) also gives rise to diffusion, 
yet smaller by a factor a;/ A, the latter can be neglected 
compared to the former in the regime we are interested 
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In order to get an efficient Monte Carlo method (im- 
portance sampling [(U))), we go over to the nonlinear ver- 
sion of (||), which keeps the trajectories ^t(z) normalized 
at all times while retaining the correct ensemble mean 
(||), see ||. Using (0), the relevant stochastic Schrodinger 
equation for Brownian motion reads 



hijj t {z) = -iH^ t (z) - i ( -m 7 A + lm{g (t)} ) q^ t {z) (8) 



+(<? - (q))zti/>t(z) 

-Re{go(t)}{(q-(q)) 2 - 
+ilm{g 1 (t)} (qp- {qp) 



((q-(q)f))Mz) 
-m(q)q - (q)p) ip t (z) 



Normalized quantum trajectories i/)t(z) whose ensem- 
ble mean gives the desired reduced density operator ac- 
cording to (||) can now be propagated using (||), where 
(q) = -faiq), a quantity which has to be determined nu- 
merically along with i/)t(z) (very often the replacement 
m(q) « (p) turns out be a good approximation). 

In (||), the time dependent coefficients quickly ap- 
proach their asymptotic values go(t) — > m 7, feT — | mj A 
and Im{gi(t)} — > — ^ for times larger than the environ- 
mental memory time. After this initial slip t ^> A -1 , (||) 
becomes 



hip t (z) 



-iH<ip t (z) 
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~lp \QP ~ (QP) + m(q)q - (q)pj ip t (z). 

We now highlight the power of our stochastic 
Schrodinger equation for Brownian motion (||) by in- 
vestigating the dynamics of a driven, damped, nonlin- 
ear, noisy system, the Duffing oscillator, where H = 
\p 2 + jQ 4 ~ \q 2 + 9Q c °s(t)i nere coupled to a heat bath 
at temperature T. This system has been studied before 
using the master equation (^) (see |ll|] and references 
therein), including a straight numerical solution which 
requires the propagation of a huge matrix. In our new 
approach, one propagates pure states i/jt(z) according to 
(||), a great reduction in resources, with the need, how- 
ever, to solve (H) many times in order to evaluate the 
mean values. For Lindblad master equations, the power 
of quantum trajectory methods for investigating classi- 
cally chaotic dissipative systems was shown in (see 
also U). 

We use parameters g = 0.3 with a damping rate 
7 = 0.25, thus the classical problem is chaotic ||l3| . The 
environment is furthermore characterized by kT = 0.3, 
and a cutoff frequency A = 5. With ft of the order 
10 -2 and smaller (see various choices of % below), the 
parameters are in the required regime. As initial con- 
dition we choose a standard coherent state located at 
(g)=0.1,(p) = 0.1. 

In Fig.l we show the ensemble mean M\W z (q,p, t = 4)] 
over 1000, 5000, and 10000 Wigner functions of pure state 
trajectories ipt(z) obtained solving (||) numerically up to 
a time t = 4. According to our construction, this quan- 
tity converges to the Wigner function of the reduced den- 
sity operator for many realizations. Here we have cho- 
sen ft = 0.01, a phase space area corresponding approx- 
imately to the extension of the wave packets shown in 
Fie.2. 
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FIG. 1. Contour plots of the Wigner function W(q,p, t = 
4) of the reduced density operator of the thermal Duffing oscil- 
lator with h = 0.01 (for the phase space area corresponding to 
this h see Fig. 2). The contour plots show the ensemble mean 
over 1000, 5000, and 10000 Wigner functions W z (q,p,t = 4) 
of individual quantum trajectories obtained solving the quan- 
tum Brownian motion stochastic Schrodinger equation (|). 

In Fig. 2 we show contour plots of Wigner functions 
W z (q,p,t = 4) of four realizations of (^), many of which 
add up to the Wigner function of the desired reduced 
density matrix shown in Fig.l. One can see clearly that 
these individual Wigner functions are well localized in 
phase space compared to the phase space spread of the 
ensemble, even for this classically chaotic system. 
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FIG. 2. Contour plots of Wigner functions W z (q,p, t = 4) 
of four individual quantum trajectories obtained solving the 
quantum Brownian motion stochastic Schrodinger equation 
(^) for the thermal Duffing oscillator. Individual trajectories 
remain well localized in phase space with respect to the overall 
spread of the ensemble mean, even for this classically chaotic 
system. The chosen value of h = 0.01 is slightly smaller than 
the phase space area covered by these states. 

This remarkable feature of the quantum Brownian mo- 
tion stochastic Schrodinger equation (||) is highlighted 
again in Fig. 3, where we show the mean position spread, 
M[Aq] = M[y ((q — {q)) 2 )] and the mean uncertainty 
product M[AqAp/h] in units of Ti of individual trajec- 
tories as a function of time for three different choices of 

h. 
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FIG. 3. Localization property of the QBM stochastic 
Schrodinger equation. Individual runs are well localized in 
phase space, the localization being stronger the smaller %: 
(a) the average position spread M[Aq] — M[^/ {(q — {q)) 2 )] 
of solutions of the QBM stochastic Schrodinger equation for 
the choices h = 0.01 (solid line), Ti — 0.005 (dashed line), and 
ft = 0.001 (dotted line). Fig. (c) shows the mean uncertainty 
product M[AqAp]/h, which remains of the order one almost 
independently of %. Thus the quantum trajectories remain 
almost minimum uncertainty wave packets for all times. 

The quantities shown in Fig. 3 can only be given sense 
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in the framework of quantum trajectories, they have no 
meaning from a density operator point of view as they 
are the ensemble mean over an expression non-quadratic 
in ijjt{z)- It is apparent from Fig. 3 that individual tra- 
jectories are well localized in phase space for all times, 
the localization being stronger the smaller Ti. As can be 
seen, our quantum trajectories remain almost 'classical' 
states, yet recover the fully quantum master equation (||). 
Thus, the representation (||) expresses the reduced den- 
sity operator of quantum Brownian motion explicitly as 
a mixture of almost 'classical' states. 

The observed localization property of QSD is well 
known in the Markov case and has been studied for in- 
stance in [0 . Here we see that similar properties hold for 
the generalized non-Markovian QSD equation (JsJ) which 
has now been applied to quantum dynamics beyond the 
class of Lindblad master equations. As in the Markov 
case, the localization property can be exploited to fur- 
ther reduce the numerical effort. 

Finally, let us briefly address the connection between 
our approach and the widely used QBM master equa- 
tion (||). Since a quantum trajectory approach strictly 
preserves positivity of the reduced density operator, our 
QBM stochastic Schrodinger equation (||) cannot be iden- 
tical to (||) in the mean, as the latter is known to vio- 
late positivity on short time scales. Taking the ensemble 
mean M [. . .] in (||) with (||) analytically, we were able to 
show in H that in the regime considered in this Letter, 
the evolution of the ensemble mean ([|) is well described 
by the master equation 

Tip = -i[H, p] - i Qm 7 A + lm{0o (*)}) [q\ p] (10) 
+ilm{#i (t)}[q, {p, p}] - Re{g (t)}[q, [q, p]] , 

which reduces to (||) for times larger than the environ- 
mental memory time, t ^> A -1 due to the asymptotics of 
the coefficients go(t),gi(t). Thus, apart from an initial 
slip on the environmental memory time scale A -1 , our 
approach recovers (^|) in the mean. It is known in the 
case of the exact master equation for a damped harmonic 
oscillator jl5| that such time dependent coefficients may 
ensure the positivity of the reduced density operator for 
non-Lindblad master equations, a result that is here sup- 
ported for general system Hamiltonian H(q,p). 

To conclude, we have presented the stochastic 
Schrodinger equation for Brownian motion. It is compat- 
ible with the standard QBM master equation yet allows 
to compute states rather than a matrix, a huge reduc- 
tion in resources, which becomes even more relevant for 
QBM in more than one space dimension. Individual tra- 
jectories are well localized in phase space, the localization 
being stronger the smaller %. Thus, in (H), the reduced 
density operator for Brownian motion is explicitly repre- 
sented as an ensemble of almost 'classical' states. 
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